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ABSTRACT 

In a universe reionized in patches, the Doppler effect from Thomson scattering off 
free electrons generates secondary cosmic microwave background (CMB) anisotropies. 
For a simple model with small patches and late reionization, we analytically calculate the 
anisotropy power spectrum. Patchy reionization can, in principle, be the main source of 
anisotropies on arcminute scales. On larger angular scales, its contribution to the CMB 
power spectrum is a small fraction of the primary signal and is only barely detectable 
in the power spectrum with even an ideal, i.e. cosmic variance limited, experiment and 
an extreme model of reionization. Consequently patchy reionization is unlikely to affect 
cosmo logical parameter estimation from the acoustic peaks in the CMB. Its detection 
on small angles would help determine the ionization history of the universe in particular 
the typical size of the ionized region and the duration of the reionization process. 

Subject headings: cosmic microwave background - cosmology: theory - intergalactic 
medium - large-scale structure of universe 



1 Alfred P. Sloan fellow 



1 



1. Introduction 

It is widely believed that the cosmic microwave back- 
ground (CMB) will become the premier laboratory for 
the study of the early universe and classical cosmol- 
ogy. This belief relies on the high precision of the 
upcoming MAPQ and Planck Surveyor^ satellite mis- 
sions and the high accuracy of theoretical predictions 
of CMB anisotropies given a definite model for struc- 



ture formation ( Hu ct al 1995). To realize the potential 



of the CMB, aspects Of structure formation affecting' 
anisotropies at only the percent level in power must be 
taken into account. 

A great uncertainty in models for structure formation 
is the extent and nature of reionization. Fortunately, 
this uncertainty is largely not reflected in the CMB 
anisotropies due to the low optical depth to Thomson 
scattering at the low redshifts in question. Reioniza- 
tion is known to be essentially complete by z ~ 5 from 
the absence of the Gunn-Pctcrson effect in quasar ab- 
sorption spectra ( Gunn fc Petcrson| 1965). Significant 



reionization before z ~ 50 will be ruled out once the ten- 
tative detections of the CMB acoustic peaks at present 



are confirmed (Scott, Silk & White 1995). It should be 
possible to deduce the reionization redshift through 
CMB polarization measurements ( ^ialdarriaga, Spergel 
k Seljak| 1997). 

Nevertheless, the duration of time spent in a par- 



tially ionized state will remain uncertain. Moreover 



as emphasized by Kaiser (1984) secondary anisotropies 



generated by the Doppler effect in linear perturbation 
theory are suppressed on small scales for geometric rea- 
sons (gravitational instability generates potential flows, 
leading to cancellations between positive and negative 

affects which are not 



Do ppler shifts). — Higher order enects wnicn are 
gen erally included in the theoretical modeling of CMB 

anisotropies are likely to be the main source of sec- 
ondary anisotropies from reionization below the degree 
scale. Such effects rely on modulating the Doppler ef- 
fect with spatial variations in the optical depth. Incar- 
nations of this general mechanism include the Vishniac 
effect from linear density variations (Vishniac 1987), the 
kinetic Sunyaev-Zel'dovich effect from clusters ( {Bunyaev 
k, Zel'dovich 1970), and the effect considered here: the 
spatial variation of the ionization fraction. 

Reionization commences when the first baryonic ob- 
jects form stars or quasars that convert a part of the 
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nuclear or gravitational energy into UV photons. Each 
such source then blows out an ionization sphere around 
it. Before these regions overlap is a period when the 
universe is ionized in patches. The extent of this period 
and the time evolution of the size and number density of 
these patches depend on the nature of the ionizing en- 
gines in the first baryonic objects. Theories of reioniza- 
tion do not give robust constraints (c/. Tegmark 1994; 

;es 1996; Aghanim et al. 1996; Loeb 1997; Haiman & 
Loet 1997; [Silk fc Rees| 1998; [Haiman fc Loeb| 1998). 



We therefore take a phenomenological approach to 
studying the effects of patchy reionization on the CMB. 
We introduce a simple but illustrative three parameter 
model for the reionization process based on the redshift 
of its onset Zj, the duration before completion 5z, and 
the typical comoving size of the patches R. It is then 
straightforward to calculate the CMB anisotropies gen- 
erated by the patchiness of the ionization degree of the 
intergalactic medium. 

We find that only the most extreme models of reion- 
ization can produce degree scale anisotropies that are 
observable in the power spectrum given the cosmic vari- 
ance limitations. A large signal on degree scales requires 
early ionization, Zi J> 30, long duration, 8z ~ z%, and 
ionization in very large patches, R ;> 30Mpc. Thus the 
patchiness of reionization is unlikely to affect cosmo- 
logical parameter estimation from the acoustic peaks in 



the CMB (Jungman et al. 1996; Zaldarriaga, Spergel & 
Seljakj 1997; [Bond, Efstathiou fc Tegmark] 1997). 



On the other hand, the patchy reionization signal on 
the sub-arcminute scale can, in principle, surpass both 
the primary and the secondary Vishniac signals. These 
may be detectable by the Planck Surveyor and up- 



coming radio interferometry measurements (Partridge 



et al.| 1997) if point sources can be removed at the 
AT/T ~ 10~ 6 level. 

An explicit expression for the CMB anisotropies 
power spectrum generated in a universe reionized in 
patches is given § [2.l| . Simple order of magnitude esti- 
mate of the anisotropy from patches is given in §|2.2|. In 



{2.2 we give a rigorous definition of our three-parameter 
reionization model, and calculate the patchy part of the 
power spectrum. We discuss illustrative examples in 



§2.4 



2. CMB power spectrum 
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2.1. Explicit expression 



2.2. Order of magnitude estimates 



Temperature perturbations A = ST/T are generated 
by Doppler shifts from Thomson scattering. For small 
optical depths 



A = — / d\ ■ vaxnxe- 



(1) 



All quantities here are in physical units. The integral 
is along the line of site, v is the peculiar velocity of 
matter, c = 1, a? is the Thomson cross section, n is the 
number density of free and bound electrons, x e is the 
local ionization fraction. 

To evaluate equation (Q) explicitly one must specify 
the cosmological model. For simplicity, we take a uni- 
verse with critical density in matter throughout; we de- 
scribe the generalization to an arbitrary FRW universe 
in § |2.3| . We furthermore use comoving coordinates x 
and conformal time 77 = f(l + z)dt = (1 + z)~ 1/,2 77o, 
where 770 = 2/Hq is the present particle horizon. We 
observe at x = and conformal time 770 along the di- 
rection of a unit vector 7; light propagation is given by 
x = 7(770 — 77). Equation (^) can be written as 

-^7 • v[tj, 7(770 - 7i)]x e [77, 7(770 - 77)]. 

(2) 

Here To = ctttio^o is the optical depth to Thomson scat- 
tering across the present particle horizon. 

The scales contributing to the peculiar velocity are 
still in the linear regime, therefore 



A(7) 



-70770 



v ( 7 ?> x ) = — v ( x )> 



(3) 



where v(x) is the peculiar velocity today. The final ex- 
plicit expression for the CMB temperature perturbation 
generated during reionization is 

A (7) = -n^o J -^i- v [i(vo-v)}xe[nn(va-v)}- (4) 

The correlation function of the temperature perturba- 
tions is defined as 



C{6) = (A( 7 i)A( 72 )) ^ 



(5) 



7l -72— cos 9 ■ 

With temperature perturbations given by equation (^) , 
this becomes 

J 'II J '12 

x x e (77i,xi)a; e (772,X2)), (6) 
where we denote Xi = 71(770 — 771) and x 2 = 72(770 — 772)- 



Consider the following patchy reionization scenario. 
The universe was reionized in randomly distributed 
patches with a characteristic comoving size R. The 
patches appeared at random in space and time. Once 
a reionized patch appears, it moves with matter. The 
average ionization fraction, that is the filling fraction of 
fully ionized patches, grows monotonically from X e = 
at high rcdshifts to X e — 1 at low redshifts. We consider 
late reionization (optical depth to Thomson scattering 
is small) and small patches (smaller than the charac- 
teristic length scale of the peculiar velocity field). We 
assume that reionization occurred at redshift Zi, and 
the patchy phase duration is given by 5z. 

The angular scale of the patchy signal is given by 
the ratio of the size of patches to the distance to them 
in comoving coordinates, 9 ~ R/(f]o — Assuming 
that the patches are uncorrelated, the spectrum of fluc- 
tuations should be white noise above this scale which 
agrees with the exact result as we shall see (eq. [p"8|). 

The rms CMB temperature fluctuation A on scales 
9 due to the patchiness can be estimated as follows. 
Since, by assumption, different patches are indepen- 
dent, A ~ A fl / 2 A p . Here N is the number of patches 
on a line of site, A p is a temperature fluctuation from 
one patch, A p ~ T p v{zi). Here v(z) = (1 + z) _1 / 2 t;(0) 
is the rms peculiar velocity at redshift z, and t p is the 
optical depth for one patch, t p ~ (1 + Z{) 2 axnaR. The 
number of patches N ~ Stj/R ~ (1 + Zi)^ 3 ^ 2 Szrjo/R. 
Collecting all the factors, we get the following estimate 
for the rms anisotropies from patches 

A ~ ro <t, 2 ) 1/2 (i?/77o) 1/2 (l + Zi) 3/4 (Szy/ 2 , (7) 

which again agrees with the exact result (eq. , up to 
a dimensionless multiplier). 

2.3. Power spectrum 

We can factor the general expression for the temper- 
ature correlation (0) as 



C(fi) 



2 4 
T oVo 



dm f dr] 2 / A . v . 
-ZT / — 3-(7i-v(xi) 72 -v(x 2 )) 



'/1 



X (x e (77l,Xi)x e (772,X2)) . 



(8) 



This assumes that x e and v are independent random 
fields. This is not strictly correct. The ionization frac- 
tion x e must be determined by the density perturba- 
tion 5, and the density perturbation is not indepen- 
dent of the peculiar velocity (for example in the linear 
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regime S — — i?yV • v). However, the ionizing radia- 
tion is presumably coming from strongly nonlinear ob- 
jects, where first stars or quasars are lightening up. At 
high z, the length scales where the density is nonlinear 
are -C lOMpc comoving, which is much smaller than 
the length scales contributing to the peculiar velocity. 
Under the assumption of scale separation, velocity and 
density (and hence x e ) are indeed independent. 

The correlation function for the local ionization frac- 
tion (x e x e ) is not known. Our model parameterizes the 
correlation function through the the patch size R and 
a mean (cosmic time - dependent) ionization fraction 



(a; e (r/i,Xi)x e (?72,X2)) = X e (r)i)X e (r) 2 ) 
+ [X e (r) mil 



X e { m )X e {r l2 )}e- i ^ 1 - 



0) 



Here r/ m j n = min(?7i, 772). The Gaussian function is cho- 
sen for simplicity; it could have been any function of 
the separation x which equals 1 at x — and gradually 
turns to zero at x > R. For the mean ionization frac- 
tion X e we assume a change from to 1 at a redshift 
Zi, with the transition occurring in a redshift interval 
Sz. We also assume Sz <C z, (this is true in all of the 
models of reionization that we are aware of). 

CMB anisotropies generated (and erased) due to the 
spatially constant part of the correlation function (^) 
are obviously the same as in the model with a uniform 
time-dependent reionization. The spatially varying part 
is responsible for generating new anisotropies; its con- 
tribution to erasing the primary anisotropies is negligi- 
ble. The anisotropy suppression is mainly determined 
by the total optical depth to Thomson scattering and is 
insensitive to the small-scale structure of the ionization 
fraction x e (j], x). 

The CMB correlation function due to the patchy part 
only is 



C {p) {9) 



2 4 
T oVo 



no 



drji 



dm 
V2 



(10) 



x (71 • v(xi)72 • v(x 2 )) e 



(*i-*2r 

2R 2 



where we denote I12 = X e (r) m i n ) — X e (rji)X e (r]2)- The 
correlation function is non-negligible only for |xi — x 2 | <; 
R. By assumption, R is much smaller then the char- 
acteristic scale of the peculiar velocity field. Also 
|x! — x 2 | <J R requires that the lines of sight 71 and 
72 be nearly parallel. Then 



(71 • v(xi)72 • v(x 2 )) 



(11) 



where (n 2 ) is the mean squared peculiar velocity today. 
For Zi ^> 1, the integral ([l0]) is dominated by small 
conformal times, and we have |xi— x 2 | 2 « 9 2 (i] — rji) 2 + 
{Vi - V2) 2 - Then 



dm f°° dm T (12-71 ) 2 , . 
< I -S^e -^-.(12) 

'0 '/1 Jo 

We assume that rji6z/(l + Zi) ^> R (with m = rj(zi)) 
and that during the patchy phase, z; > z > zi — Sz, 
the ionization fraction X e grows linearly from to 1 
such that eventually both hydrogen and helium arc fully 
ionized. Then 



C (p) (<9) =Ae 
where the characteristic angular scale is 
R R (1 + z,) 1 / 2 



Vo-Vi Vo (1 + Zi) 1/2 - 1' 
and the amplitude is 



.4 



R 



• 2 (11 2 ) —Sz(l + Zi) 



3/2 



(13) 



(14) 



(15) 



36 " 1 ' rjo 

Note that a critical matter-dominated universe is as- 
sumed in this expression. To generalize this result re- 
place 7/0 — r\i by the comoving angular diameter distance 
in equation ( |l4| ) and a factor of (1 + Zi) in equation (15) 
with the appropriate velocity growth factor. 

The power spectrum is given by the spherical har- 
monics decomposition 



(p) 



2?r / dcos6Pi(cos6)w p (9). 



For I > 1, 

C ; (p) w 2tt / 6d6J {W)w p {9) = 2irA6 2 l e-°- Si ir 



The power per octave is 
2r»(p) 



l 2 C, 



2- 



Al 



2a2, 



(16) 



(17) 



(18) 



The anisotropy power reaches the maximal value 



,i 2 c[ v) . 

1 2tt ' 



' 27r T g(v*)*6z(l + Zif/2, (19) 



at 



18e 
R 



[l-(l + z t )-^ 2 ] 



(20) 
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2.4. Discussion 

The signal from patchy reionization in our model 
depends on four quantities: the rms peculiar velocity 

1/2 

^u 2 ) today, the redshift of reionization zi, its du- 
ration 5z and the characteristic comoving size of the 
patches R. The structure formation model specifies the 
power spectrum of fluctuations which in turn tells us the 
rms peculiar velocity. Let us now consider the patchy 
reionized signal in the context of a specific model for 
structure formation. 

For illustrative purposes, let us consider a cold dark 
matter model with h = 0.5, Ofe = 0.1, and a scale- 
invariant n = 1 spectrum of initial fluctuations. Nor- 
malizing the spectr um to the COB E detection via the 
fitting formulae of Bunn fc White (1997) (their equa- 



tions [17]- [20]) and employing the analytic fit to the 
transfer function of Eiscnstcin fc Hu| (1998) (their equa- 
tions [15]-[24]) we find an rms velocity of (u 2 ) 1 ^" = 
3.9 x 10~ 3 . With the present optical depth of tq = 
0.122£!b/i = 0.0061, we have a maximal anisotropy of 

(^W = 2.41 x 10- 15 ^Uz(l + z t fl\ (21) 



at 



16958 
i?/Mpc 



Mpc 



[1- (1 + ^-1/2]. 



(22) 



The power spectrum of the model in principle also 
tells us the remaining parameters of the ionization: its 
redshift Zi, duration 5z and typical patch size R. Un- 
fortunately, these quantities depend on details of the 
cooling and fragmenting of the first baryonic objects 
to form the ionizing engines. We therefore consider 
5 <J Zi 50 which spans the range of estimates in 



the liter ature ( Tegmark 1994; Rees 1996; Aghanim ct 
19 96; |Loeb| 1997; |Haiman fe Loeb| 1997; |Silk fc Ree; 
1998; Haiman fc Loeb 1998). Reionization, once it com 



mences, is generally completed in a time short compared 
with the expansion time at that epoch 8z/(l+Zi) < 1 by 
the coalescence of patches that are small compared with 
the horizon at the time Rjr\i <C 1 at the time. Again 
the exact relations depend on the efficiency with which 
the first objects form and create ionizing radiation (see 
e.g. fTegmark| 1994). 

Let us consider an extreme example of z, ; = 10, Sz = 
3, R — 20Mpc. Then the maximal power is « 5.3 x 
10 -12 at I ss 590, the primary signal at these scales is w 
3 x 10 -10 - the contribution of the patchy reionization is 
small in comparison (see Fig. [j]) . However, in light of the 
high precision measurements expected from the MAP 
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Fig. 1. — CMB anisotropy power spectra in a CDM 
model with extreme patchiness. Shown here are the pri- 
mary anisotropy and the patchy reionization anisotropy, 
eq. (£§) with Zi = 10, 5z = 3, R = 20Mpc. These sig- 
nals are compared with the cosmic variance of the pri- 
mary anisotropy and the noise of the MAP satellite (in 
logarithmic bins). 

and Planck satellites such a signal is not necessarily 
negligible. The ultimate limit of detectability through 
power spectrum measurements is provided by so-called 
cosmic variance. This arises since we can only measure 
2£+l realizations of any given multipole such that power 
spectrum estimates will vary by 



SC,= 



21 + 1 



-a 



(primary) 



(23) 



Detection of a broad feature such as that from patchy 
reionization is assisted in that we may reduce the cos- 
mic variance by averaging over many ^'s. We show an 
example of this averaging in Fig. [j] (lower left boxes). 
In this model, the patchy reionization signal can be de- 
tected at the several a level if cosmic variance were the 
main source of uncertainty. Of course a realistic ex- 
periment also has noise and systematic errors. We also 
show the noise error contributions expected from the 
MAP experiment in Fig. [I]. 

An important additional source of uncertainty is pro- 
vided by other unknown aspects of the model. In- 
deed it is hoped that the CMB power spectrum can be 
used to measure fundamental cosmological parameters 
to high precision. Excess variance from patchy reion- 
ization can in principle cause problems for cosmological 
parameter estimation from the CMB if not included in 
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the model. It would remain undetected and produce 
parameter misestimates if its signal can be accurately 
mimicked by variations in the other parameters. For- 
tunately, the angular signature we find here - £ 2 white 
noise until some cut off due to the patch size - does not 
resemble the signature of other cosmological parameters 
which alter the positions and amplitudes of the acoustic 



peaks (see Bond, Efstathiou & Tegmark 1997; Zaldar 
ria; ;a, Spergel fc Seljak] 1997). Coupled with the small 



amplitude of the effect on the 10 arcminute to degree 
scale for even this extreme model, it is unlikely that 
patchy reionization will significantly affect parameter 
estimation through the CMB. 

We have called the (z^ = 10, Sz = 3, R = 20) model 
extreme, because of the size of patches; the reionization 
redshift and duration would be considered reasonable 
by a number of theories. For example the early quasar 
model of Haiman & Loeb (1998) does predict Zi ~ 10 
and 8z ~ 3. However, their "medium quasar" emits 
only ~ 10 67 ionizing photons during its life time. These 
photons cannot ionize a bubble larger than R ~ IMpc 
comoving. 

Perhaps more interesting is the case where reioniza- 
tion takes place at a higher redshift with for example 
Zi = 30, Sz = 5, R = 3Mpc. The reduction in the 
patch size causes the signature to move to smaller an- 
gles where the primary signal is negligible due to dissi- 
pational effects at recombination. The increase in the 
optical depth at this higher redshift is counterbalanced 
by the reduction in the rms fluctuation due to the num- 
ber of patches along the line of site such that the am- 
plitude of the signal increases only moderately. Here 
the maximal power is w 6.2 x 10~ 12 at I w 4650 (see 
Fig. ||) . Patchy reionization effects exceed the Vishniac 
signal at these scales (~ 3 x 10~ 12 ) which is believed 
to be the leading other source of secondary anisotropics 



(Hu fc White 1996) 



Although the morphology and amplitude of the patchy 
reionization and Vishniac signals are similar, the Vish- 
niac effect is fully specified by the ionization redshift 
and the spectrum of initial fluctuations and hence may 
be removed once these are determined from parame- 
ter estimation at larger angular scales. Likewise, since 
the rms peculiar velocity (v) and the ionization redshift 
Zi will be specified by the large scale observations, the 
amplitude of the signal can be used to estimate the du- 
ration of reionization Sz and its angular location the 
typical comoving size of the bubbles R. 

In summary, the patchiness of reionization leaves a 
potentially observable imprint on the CMB power spec- 




Fig. 2. — CMB anisotropy power spectra in a CDM 
model with early reionization. Shown here are the pri- 
mary anisotropy suppressed by rescattering and the 
patchy reionization anisotropy, eq. ( |l8| ) with zi = 30, 
Sz = 5, R = 3Mpc. These signals are compared with 
the cosmic variance of the primary anisotropy achiev- 
able by an ideal experiment in the absence of galactic 
and extragalactic foregrounds. 



trum, but one that is unlikely to affect cosmological 
parameter estimation from the acoustic peaks in the 
CMB. We show how the signature scales with the gross 
properties of reionization - its redshift, duration, and 
typical patch size. Observational detection of this sig- 
nature would provide useful constraints on the presently 
highly uncertain reionization scenarios but will likely 
require experiments with angular resolution of an ar- 
cminute or better and foreground subtraction at better 
than the 8T/T ~ 10~ 6 level. 
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